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THE  EFFECT  OF  WEAK  VISCOUS  INTERACTION 


ON  THE  RESISTANCE  OF  A WING  SECTION 


V.  Ya.  Ivanov  and  V.  M.  Kovalenko 


As  is  known  the  lines  of  flow  in  the  boundary  layer  separate 
from  the  streamlined  surface  as  a result  of  the  light  braking  action 
on  the  part  of  the  wall,  which  leads  to  a change  in  pressure  dis- 
tribution along  this  surface.  In  other  words,  the  boundary  layer 
governed  by  the  external  flow  exerts  an  inverse  effect  on  this  flow 


- what  we  have  is  the  so-called  effect  of  viscous  interaction. 


From  the  practical  standpoint  the  estimation  of  the  effect  of 
the  interaction  of  the  boundary  layer  with  the  external  flow  on  the 
actual  wing  section  is  of  considerable  interest.  Such  an  inter- 
action is  manifested  in  the  change  in  the  wave  drag  of  the  profile 
and  in  its  friction  drag.  The  actual  pressure  distribution  along 
the  wing  section  surface  with  a two-dimensional  motion  of  the 
viscous  gas  flow  coincides  with  the  pressure  distribution  during 
a smooth  flow  of  an  ideal  gas  around  a body,  formed  by  the  builtup 
of  the  displacement  thickness,  calculated  according  to  the  actual 
pressure  distribution  [1],  of  the  boundary  layer  on  the  wing  sec- 

tion and  on  both  sides  of  the  zero  line  of  flow  in  its  wake.  With 
this  approach  the  problem  of  determining  the  actual  pressure  dis- 
tribution is  solved  using  the  method  of  successive  approximations. 
The  numerical  calculations  have  shown  that  even  the  first  approxima- 
tion enables  us  to  obtain  the  principal  quantitative  portion  of 
the  correction  which  takes  into  account  the  inverse  effect  of  the 


i 


boundary  layer  on  the  theoretical  distribution  of  pressure. 

The  effect  of  viscous  interaction  at  hypersonic  speeds  has  been 
studied  extensively  in  [2,  3 ] i less  attention  nas  been  devoted  to 
the  problem  dealing  with  the  effect  of  a weak  interaction  on  the 
drag  of  a wing  section  at  moderate  supersonic  speeds,  in  particular 
with  the  presence  of  a mixed  boundary  layer  on  such  a profile. 

This  question  has  been  formulated  most  extensively  in  [4,  5]- 
We  will  use  the  following  designations  in  this  work:  x,  s - 
coordinates  along  the  chord  and  along  the  profile  generatrix,  res- 
pectively, x=x/b;  b - length  of  the  profile  chord;  c - profile 

thickness,  c=c/b;  l - extent  of  the  transition  zone  of  the 

t , n 

boundary  layer  from  the  laminar  state  into  turbulent  state,  l = 

t 9 n 

Z T n /b ; 6j,  - displacement  thickness  and  thickness  of  the 

impulse  loss  of  the  boundary  layer;  p,  y,  and  T - density,  viscosi 
ty  factor,  and  gas  temperature; 


Re*  = 


i 


. Rex  = 


“ft*  Pa 
P« 


Reynolds  numbers, 

_ T 

M - Mach  number;  T^=  ^ - temperature  factor;  c^,  - local  coef- 
ficient of  friction  drag;  cx  - coefficient  of  profile  resistance. 
The  subscripts:  6 pertains  to  the  parameters  on  the  external  line 
of  the  boundary  layer;  w - to  the  parameters  on  the  wall;  e - to 
the  parameters  on  the  heat-insulated  wall;  » - to  the  parameters 
of  the  incident  flow;  n,  t - laminar  and  turbulent  layers,  respec- 
tively; and  n/i.  - plate. 

In  this  work  we  Investigate  the  effect  of  viscous  interaction  on 
a symmetrical  parabolic  profile  (Pig.  1).  It  is  assumed  that  the 
shock  wave  arising  on  a sharp  leading  edge  of  the  profile  is  recti- 
linear. This  fact  excludes  from  the  examination  the  cases  of 
relatively  thick  profiles  at  low  supersonic  speeds.  It  is  also 
presupposed  that,  in  contrast  to  works  [4,  5],  the  transition  of 
the  boundary  layer  from  laminar  to  turbulent  state  occurs  not  as 
a shock  wave  but  occupies  a certain  region.  Furthermore,  the 
turbulent  boundary  layer  Is  calculated  by  a more  perfect  method 
and  the  estimation  of  the  additional  wave  drag  of  the  friction  was 
accomplished  over  a wide  range  of  change  in  the  basic  parameters 


* • • ' ■ s r ' 1 r 
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of  the  incident  flow  and  the  geometrical  dimensions  of  the  profile. 

The  distribution  of  speeds  along  the  initial  profile  can  be 
written  in  the  form  of  the  functional  relationship 


«4o  = / IP  (I).  P (0)J.  (1) 

Here  um^um/uo,:  uto  - speed  at  the  outer  limit  of  the  boundary 
layer  without  the  consideration  of  the  interaction;  u«,  - velocity 
of  the  incident  flow;  S(x)  - angle  between  the  direction  of  the 
incident  flow  and  the  tangent  to  the  initial  contour  of  the  profile 
at  its  every  point;  B(0)  - initial  rotation  angle  of  the  flow. 

The  velocity  of  the  external  flow  u^  + disturbed  by  the  boundary 
layer,  in  accordence  with  what  was  said  above,  can  be  written  in 
the  form 

«fr+  = / IP W + Ap,  P (0)  + Ap0],  (2) 

where  due  to  the  smallness  of  AB 

Ap*tgAp  = ^.  (3) 


The  assumption  regarding  the  smallness  of  AB , which  is  valid 
for  the  points  sufficiently  removed  from  the  leading  edge,  is  a 
necessary  condition  of  the  weak  interaction  j . However, 

close  to  the  leading  edge  the  derivative  d6 j/ds  is  large  and  there 

is  strong  interaction  . But  since  the  area  of  strong 

interaction  is  very  small  on  thin  tapered  bodies,  its  contribution 
to  the  overall  effect  is  negligible.  In  view  of  this  fact,  when 
calculating  the  region  close  to  the  leading  edge,  as  was  done  in 
[5]»  we  used  the  condition 


dfi, 

d*  0 <7  <0.05 


dfi, 

ds  7 0,05 


— Ap0, 


(*0 


which  allows  us  to  calculate  this  area  in  the  approximation  of  a 
weak  interaction. 

The  distribution  of  the  velocities  and  pressures  along  the  profile 
contour  was  calculated  according  to  the  A.  Ye.  Donov's  theory  of 
third  approximation  [6],  i.  e.,  it  was  assumed  that 
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««o  = 1 + M (x)  -f  6..p*  (x)  + b'i  Ps  ( X ) + (&,  — bi)  P5  (0),  ^ 5 ) 

Pm  = a,P  (x)  + a,pV)  + a, p* > x)  + ald pJ  (0) . ( 6 ) 

In  a particular  case,  when  p(0)^0  , the  coefficients  b^b^  anc3 

a iu=0  . In  formulas  (5)  and  (6)  the  pressure  coefficient 

— Pfith  — P 

p«o=  — a — ~ » values  poo  and  qoo  are  pressure  and  velocity  head  of 

’ 00 

the  incident  flow,  respectively,  b.^  and  a ^ are  functions  of  the 
Moo  number  [6]. 

If  expand  (2)  in  the  Taylor  series,  using  relationships  (4)  and 
(5),  and  keep  the  terms  only  to  the  third  degree  of  smallness  in- 
clusively (in  this  case  ^ - value  with  the  second  degree  of  small- 
ness), equation  (2)  will  assume  the  form 

«^+  = ^+(61  + 261p(J)lgi  + Aufti  (7) 

where  value  Au^  is  determined  by  the  initial  (i.  e.,  when  x=0) 
values  of  the  angles  3(0)  and  ABo  and  is  necessary  for  a correct 
recalculation  of  the  shock  wave  with  the  consideration  of  the  inter- 
action . 

The  formula  for  the  pressure  coefficient  p^+  can  be  obtained  in 
the  same  way  with  the  consideration  of  the  interaction  of  the 
boundary  layer  with  the  external  flow: 

^6+  = P<*  + dP  + *P’  ( 8 ^ 

where 

dp  = [ aj  + 2a*p  (x)  + 3a,p*  (x)j^,  ( 9 ^ 

Pgo  is  determined  by  relationship  (6),  and  Ap  in  the  same  way  as 
the  value  Au.  by  the  values  3(0)  and  ^1 

0 ds  £=0,05. 

First,  formulas  (5)  and  (6)  are  used  to  calculate  the  external 
flow.  The  determined  distributions  of  pressure  and  velocity  along 
the  contour  allowed  us  to  determine  the  values  of  the  displacement 
thickness  of  the  boundary  layer  in  the  first  approximation.  Then, 
using  formulas  (7)  and  (8),  the  external  flow  is  calculated  in  the 
second  approximation,  taking  into  account  the  effect  of  interaction. 
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According  to  the  data  in  [4,  5],  the  estimation  of  this  effect  can 
be  limited  by  the  second  approximation. 

The  calculation  of  the  wave  drag  c for  one  side  of  the  profile, 
and  also  of  the  additional  resistance  Ac  caused  by  the  effect  of 
interaction,  was  accomplished  in  a rapid  coordinate  system  using 
the  simplified  formulas  obtained  from  the  exact  expressions  for  the 
c and  Ac  [6]  with  the  expansion  of  the  trigonometric  functions, 
which  enter  into  them,  in  a series  of  powers  of  the  corresponding 
angles  with  an  accuracy  to  the  terms  which  are  proportional  to  the 
cubes  of  the  angles  (which  is  possible  for  a slender  profile  at 
relatively  small  angles  of  attack).  These  formulas  have  the  form 


Here 


Cx,  = f + pfl  (X)  dx, 

i 

&CXU  = j*  (dp  + A p)  Pa  (x)  dx, 

b 

(x) = r>  fx)  1 + ? ± aP  w]  + if>s  w 


(+  for  the  upper  and  lower  surfaces  of  the  profile,  respectively); 
a - angle  of  attack. 

I-n  a general  case  the  estimation  of  the  effect  of  viscous 
interaction  was  accomplished  in  the  presence  of  a mixed  boundary 
layer.  The  characteristics  of  the  laminar  layer  were  calculated 
according  to  the  method  in  [7],  which  is  suitable,  however,  only 
for  the  relatively  small  gradients  of  pressure.  The  method  in  [8] 
was  used  to  calculate  the  turbulent  boundary  layer;  this  method 
can  be  used  for  a continuous  flow  around  a wing  section  with  at- 
tached shock  wave.  The  momentum  equation  is  used  to  determine  the 
local  coefficient  of  friction  drag  which  can  be  written  in  the 
form 

clt = — I iniL* ^ 4.  r d Ke “ j.  1^1 

T'  T*-  “*  !■  (13) 


where 


"a  f'4  6 1 


= Rew  ( clt , Mft,  Tw). 
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"^■^^UKUHP 


m — wm  u 


Thickness  of  the  momentum  loss  6^  is  determined  using  the  data 
for  a plate  and  the  form-parameter  H is  determined  according  to  the 
formula 


H = 1.285  Tu  + 0,2  U5,r„). 


(15) 


As  the  boundary  condition  we  used  the  condition  where  in  the 
first  calculation  point,  which  5s  relatively  close  the  leading 
edge  (when  x=0.03), 


cit\x-o,o:i  . (M*  |t=o,03.  tfe,  |;«0.o3,  TJ. 


(16) 


In  the  case  of  a mixed  boundary  layer  the  determination  of  the 
boundary  condition  depends  on  the  assumptions  regarding  the  relative 
length  of  the  transition  zone  1 T n<  If  we  assume  that  the  transi- 
tion occurs  instantly  (i.  e.,  1^  n=0),  the  local  coefficient  of 
resistance  of  the  turbulent  friction  at  the  point  of  transition  can 
be  determined  from  the  condition  of  equalitv  of  the  momentum  loss 
thicknesses  of  the  laminar  and  turbulent  boundary  layers.  In  the 
case  when  the  length  of  the  transition  zone  has  a finite  dimension 
(.1  >0),  the  boundary  condition  for  equation  (13)  can  be  determined 

by  introducing  certain  assumptions.  The  displacement  thickness  of 
the  boundary  layer  at  the  end  of  the  transition  zone  was  assumed 
to  be  equal  to  the  half-sum  of  the  displacement  thicknesses  of  the 
laminar  and  turbulent  boundary  layers,  which  were  calculated  for 
the  point  on  the  profile  which  is  located  in  the  middle  of  the 
transition  zone;  in  this  case,  the  thickness  Sj  at  this  point  is 
calculated  from  the  condition  ^2n=^2T  ^see  Fig.  2).  The  measure- 
ment results  obtained  for  the  displacement  thickness  on  a flat 
plate  in  the  transition  zone  can  serve  as  a certain  confirmation 
of  this  assumption  [9]. 

Equation  (13)  was  solved  numerically  according  to  the  Runge- 
Kutta  method  with  a fixed  step,  which  equals  to  0.01  of  the  chord. 
After  the  values  of  the  c^g  were  determined,  we  calculated  the 
local  coefficients  of  the  friction  drag  c-  , referred  to  the  con- 
ditions  in  an  incident  flow.  In  determining  the  average  friction 
coefficient  Cp,  i.  e.,  when  integrating  the  values  of  c^  along 
the  profile's  contour  and  projecting  the  friction  force  onto  the 
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direction  of  the  incident  flow,  it  was  assumed  that,  in  the  transi- 
tion zone,  the  local  friction  coefficient  changes  linearly. 

Change  in  the  pressure  of  an  inviscid  flow,  arising  due  to  the 
effect  of  the  boundary  layer,  x’ads  to  the  change  in  the  surface 
friction  determined  according  to  the  initial  pressure  distribution 
in  an  inviscid  flow,  thus 


Ac?  = Cf+  — Cf  o. 


(17) 


where  c^+  - average  friction  drag  coefficient  determined  with  the 
consideration  of  the  interaction;  CpQ  - without  the  consideration 
of  the  interaction. 

The  total  effect  derived  from  the  interaction  between  the  bound- 
ary layer  and  the  external  flow  can  be  expressed  by  means  of  the 
coefficient 

Acr.\  / . Ac*.  \ Acr  ( 18 ) 


Aa  l = ( 


I 4- 


Lro 


A Cp  ' 

CF0  ‘ 


1 


lfq 


The  calculations  were  carried  out  on  a computer  according  to 
the  established  procedure  over  a wide  range  of  change  in  the  basic 
parameters  of  the  flow  (2^M„^6,  107<Re.^  109,  0,6 ^ ^ 1 ) and  the 

geometry  of  the  profile  0^c^0,06)  . The  static  parameters  of  the 
incident  flow  (T„,  p„,  p»,  g.)  correspond  to  the  altitude  of  20  km. 
Figures  3 to  7 show  some  calculation  results  obtained  for  the  coef- 
ficients of  additional  resistance  Acx , and  A cF  , attributed  to 
the  friction  resistance  of  the  initial  contour.  The  comparison  with 
the  friction  drag  is  made  due  to  the  fact  that  the  source  of  both 
components  of  the  additional  resistance  is,  ultimately,  the  viscos- 
ity of  the  flow. 

From  an  analysis  of  the  calculation  material  we  can  conlude  that 
the  thicker  the  boundary  layer  and,  consequently,  the  greater  the 
displacement  thickness,  i.  e.,  the  larger  the  Moo  (Fig.  3)  and  T 
(Fig.  5)  numbers  and  the  smaller  the  Re*,  number  (Fig.  *0  the 
greater  the  additional  resistance  of  the  profile.  Furthermore,  for 
two  profiles  which  differ  in  thickness  and  which  are  in  the  ‘"ame 
flow,  greater  additional  resistance  arises  on  a thicker  profile 
(Fig.  6).  In  comparison  with  the  friction  drag,  for  example  at  the 

7 

numbers  Moo=3  and  Re  =10  , the  portion  of  this  resistance  comorises 
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and  it  rapidly  increases  with  an  increase  in  the  number. 

Of  practical  interest  is  the  question  concerning-  the  influence 
of  the  position  of  the  transition  zone  and  its  extent  on  the  com- 
ponents of  additional  resistance  of  the  profile.  The  calculations 
have  shown  that,  in  the  case  of  a mixed  boundary  layer,  there  are 
ranges  of  basic  parameters  of  a flow  at  which  the  interaction  may 
have  no  effect  on  the  wave  or  total  resistance  of  the  wing  section, 
while  in  other  cases,  it  can  lead  to  a decrease  in  the  wave  and 
total  resistances.  The  maximum  of  this  decrease  corresponds  to 
the  case  where  the  point  of  transition  is  close  to  the  middle  of 
the  profile;  its  value  depends  on  the  flow  parameters  (Fig.  7). 

The  extent  of  the  transition  region  virtually  has  no  effect  on  the 
magnitude  of  the  additional  wave  resistance  within  the  limits  of 
this  region  and  and  it  is  manifested  only  in  the  change  of  the 
specific  contribution  to  the  value  c,„  of  the  purely  laminar  and 
purely  turbulent  areas  of  the  flow. 

The  calculations  show  that  the  additional  resistance  of  the  wing 
section,  caused  by  the  effect  of  viscous  interaction,  is  significant 
for  important  cases  in  practice  and  it  must  be  taken  into  account. 

It  is  possible  to  assume  that  the  values  of  Ac*  , calculated  in  this 
manner,  remain  in  force  also  for  the  profile  which  is  in  the  flow 
at  a small  angle  of  attack. 
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KEY:  1)  Laminar  layer  2)  Tran- 
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Pip.  2.  Distribution  of  the  dl s 
placement  thickness  in  the  tran 
sition  zone  [9] 


Pip.  3.  Additional  total  drap  as 
a function  of  the  number  for 
a plate  and  a parabolic  profile 
with  the  thickness  c=0.06.  Re  = 


Pig.  A.  Additional  wave  drap 
and  Priction  drap  as  a function 
of  the  Rero  for  a parabolic 

profile  (c=0 .06).  Mro=6>  Tw=1.0. 

turbulent  layer. 
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lami nar 


laminar  layer 


C, 

AC,,/ 

"\ 

^ 



- 

Fig.  5.  Effect  of  the  temperature 
factor  on  the  magnitude  of  the 
additional  wave  drag  and  friction 
drag  of  a parabolic  profile 

(c=0.06),  =6 , Reoo=10^,  x,p=0. 


Fig.  6.  Effect  of  the  thickness 
of  a parabolic  profile  on  the 
magnitude  of  the  additional 
wave  drag,  Moo=6,  Re^lO  , 

Tw=l. 0 . 

turbulent  layer,  laminar 

layer. 


Fig.  7.  Additional  resistance  of  a plate  and 
two  parabolic  profiles  (c=0.02  and  0,06),  due 
to  the  effect  of  viscouse  interaction,  depend- 
ing on  the  position  of  the  point  of  transition. 
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